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ABSTRACT 

The present study examines self-similarity evolution of advection dominated accre- 
tion flow (ADAF) in the presence of a toroidal magnetic field. In this research, it 
was assumed that the angular momentum transport is due to viscous turbulence and 
a-prescription was used for kinematics coefficient of viscosity. The flow does not have 
a good cooling efficiency and so, a fraction of energy accretes with matter on central 
object. The effect of a toroidal magnetic field on such systems in a dynamical be- 
havior was investigated. In order to solve the integrated equations which govern the 
dynamical behavior of the accretion flow, self-similar solution was used. The solution 
provides some insights into the dynamics of quasi-spherical accretion flow and avoids 
many of the strictures of the steady self-similar solutions. The solutions show that 
the behavior of physical quantities in a dynamical ADAF are different from steady 
accretion flow and a disk with polytropic approach. The effect of the toroidal magnetic 
field is considered with additional variable f3[= p m ag/Pgas], where Pmag and Pgas 
the magnetic and gas pressure, respectively. Also to consider the effect of advection 
in these systems, the advection parameter / was introduced that stands for a fraction 
of energy that accretes by matter to the central object. The solution indicates a tran- 
sonic point in the accretion flow for all selected amounts of / and j3. Also, by adding 
strength of the magnetic field and the degree of advection, the radial-thickness of the 
disk decreased and the disk compressed. The model implies that the flow has differen- 
tial rotation and is sub-Keplerian at small radii and is super-Keplerian in large radii 
and that different result was obtained using a polytropic accretion flow. The obtained 
(3 parameter was used a function of position that increases by increasing radii. Also, 
The behavior of ADAF in a large toroidal magnetic field implies that different result 
was obtained using steady self-similar models in large magnetic field. 
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1 INTRODUCTION 

During recent years one type of accretion disks has been 
studied, in which it is assumed that the energy released 
through viscous processes in the disk may be trapped within 
the accreting gas. This kind of flow is known as advection- 
dominated accretion flow (ADAF). The basic ideas of such 
ADAF models have been developed by a number of re- 
searchers (e.g., Ichimaru 1977; Rees et al. 1982; Narayan 
& Yi 1994, 1995; Abramowicz et al. 1995; Ogilvie 1998; Ak- 
izuki & Fukue 2006; hereafter AF06). 

It is thought that accretion disks, whether in star- 
forming regions, in X-ray binaries, in cataclysmic variables, 
or in the centers of active galactic nuclei, are likely to be 
threaded by magnetic fields. Consequently, the role of mag- 
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netic fields on ADAF has been analyzed in detail by a 
number of investigators (Bisnovatyi-kogan & Lovelace 2001; 
Shadmehri 2004; AF06; Ghanbari et al. 2007, Abbassi et 
al. 2008). The existence of the toroidal magnetic fields have 
been proven in the outer regions of YSO discs (Greaves et 
al. 1997; Aitken et al. 1993; Wright et al. 1993) and in the 
Galactic center (Novak et al. 2003; Chuss et al. 2003). Thus, 
considering the accretion disks with a toroidal magnetic field 
have been studied by several authors (AF06; Begelman & 
Pringle 2007; abbassi et al. 2008; Khesali & Faghei 2008 
and references within; hereafter KF08). KF08 considered 
dynamic behavior of a polytropic accretion flow in presence 
of a toroidal magnetic field. In a dynamic approach they 
showed the radial behavior of the physical quantities were 
different with results achieved by those who considered the 
accretion flow in a steady self-similar methods (Shadmehri 
2004; AF06; Ghanbari et al 2007; Abbassi et al. 2008). For 
example, KF08 presented that ratio of the magnetic pres- 
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sure to the gas pressure is not constant and varies by radius. 
The results of KF08 were assembled on polytropic equation 
that implies the accreting gas has a good cooling efficiency, 
while results of some authors have shown that the behavior 
of physical quantities are very sensible to fraction of the en- 
ergy that traps within the accreting gas (AF06). So,in the 
present study it is intended to investigate dynamic behavior 
of an ADAF in presence of a toroidal magnetic field. The pa- 
per is organized as follows. In section 2, the general problem 
of constructing a model for quasi-spherical magnetized ad- 
vection dominated accretion flow will be defined. In section 
3, self-similar method for solving the integrated equations 
which govern the dynamic behavior of the accreting gas was 
utilized. The summary of the model will appear in section 



Here p is the density, v r the radial velocity, D. the angular 
velocity, M» the mass of the central object, p the gas pres- 
sure, B v the toroidal component of magnetic field, v the 
kinematic viscosity coefficient and it is given, as in Narayan 
& Yi (1995), by an a-model 

(6) 



pQ, K 



where Qk = (GM,/r 3 ) 1 '' 2 is the Keplerian angular veloc- 
ity. The parameters 7 and a are assumed to be constant 
and / measures the degree to which the flow is advection- 
dominated (Narayan & Yi 1994), and is assumed to be con- 
stant. 



2 BASIC EQUATIONS 

In this section, we derive the basic equations which describe 
the physics of accretion flow with a toroidal magnetic field. 
We use the spherical coordinates (r, 9, (f>) centred on the ac- 
creting object and make the following standard assumptions: 

(i) The accreting gas is a highly ionized gas with infinitive 
conductivity; 

(ii) The magnetic field has only an azimuthal component; 

(iii) The gravitational force on a fluid element is charac- 
terized by the Newtonian potential of a point mass, ^ = 
— GM*/r, with G representing the gravitational constant 
and M» standing for the mass of the central star; 

(iv) The equations written in spherical coordinates are con- 
sidered in the equatorial plane 9 — tt/2 and terms with any 
9 and ip dependence are neglected, hence all quantities will 
be expressed in terms of spherical radius r and time t; 

(v) For simplicity, the self-gravity and general relativistic 
effects have been neglected. 

Under the assumptions and the approximation of quasi- 
spherical symmetry and the ideal magnetohydrodynamics 
treatment, the dynamics of a magnetized accretion flow is 
described by the following equations: 
the continuity equation 

dp . 1 d 2 s n 
the radial force equation 
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dv r dv 
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p dr r 
the azimuthal force equation 
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the energy equation 

1 " dp ^ dp 
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and the field freezing equation 

dB v 1 d 
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3 SELF-SIMILAR SOLUTIONS 
3.1 Analysis 

Self-similar models have proved very useful in astrophysics 
because the similarity assumption reduces the complexity 
of the partial differential equations. Even greater simplifi- 
cation is achieved in the case of spherical symmetry since 
the governing equations then reduce to comparatively sim- 
ple ordinary differential equations. We introduce a similarity 
variable r\ and assume that each physical quantity is given 
by the following form: 



r = r (t)ri 
p(r,t) = po(t)R(v) 
p(r,t) = Po(t)P(v) 
v r (r,t) = v (t)V(r)) 
Q(r,t) = n (t)io(ri) 
B v (r,t) = b Q {t)B{r)). 



(7) 
(8) 
(9) 
(10) 

(11) 
(12) 



By assuming power-law time dependent of ro(t) = at", 
where n = 2/3, we find the following relations: 



r (t) = at 2 ' 3 
Po(t)/po(t) 
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(13) 
(14) 

(15) 

(16) 
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The above results imply that po(t) and bo(t) are dependent 
on timely behavior of po(t). So, for specifying time depen- 
dent of po(t), and then po(t) and bo(t), we introduce the 
mass accretion rate M 



M 



-47IT pV r - 
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Similar to equations (7)-(12) for the mass accretion rate we 
can write 



M(r,t) = M (t)rh(v)- 



(19) 



Under transformations of equations (7), (8) and (10), equa- 
tion (19) becomes 
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M(r,t) = [r 2 (t)po(t)vo(t)] x [-4irr, 2 R(rj)V(r,)] (20) 
in which implies 

M {t) = r 2 (t)po(t)v {t) (21) 

m{rj) = -4irrj 2 R{rj)V{rj). (22) 

Now, we consider a set of solutions that Mo(t) is a constant 
(KF08), thus we can write 



(23) 
(24) 



p (t) = (Mo/VGM^K 1 
that implies 

po(t) = (Moy/ GM */ a5 ) t ~ 5/3 
and 

&g(*)/87r = {M y/GM t /a 5 )t~ 5/3 . (25) 

Substituting equations (6)- (12) and (13)- (17) into the basic 
equations (l)-(6), the similarity equations are obtained as 

2r,\dR , R d t _ 2jA _ n (26) 
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To investigate existence of transonic point, the square 
of the sound velocity is introduced that subsequently can be 
expressed as 



a = P _ GIVLP 2 /3 
p an 



(31) 



Here, S = (P/R) 1 ^ 2 the sound velocity in self-similar flow, 
which is rescaled in the course of time. The Mach number 
referred to the reference frame is defined as (Fukue 1984; 
Gaffet & Fukue 1983) 

v r — vf _ V — nrj 



S 



(32) 



(33) 



H = 



where 

dr r 

VF = Tt= n l 

is the velocity of the reference frame which is moving out- 
ward as time goes by. The Mach number introduced so far, 
represents the instantaneous and local Mach number of the 
unsteady self-similar flow. We will consider transonic points 
of accretion flow in next subsection. 



In order to consider the strength of the magnetic field 
in the plasma, the /3 parameter is introduced that is ratio of 
the magnetic to the gas pressures 



/3(r,t) = ^ 



B 2 (r,t)/8n _ B 2 (rj) 



p(r,t) 



P(rj) 



(34) 



In completing this section, we also summarize the main 
results here. Solving equations (1), (10), (11), and (19) under 
transformations (12)- (15) in non- magnetically state, makes 
it clear that time behavior of physical quantities in the non- 
magnctically and the magnetically disk are the same. This 
result is one of the strictures of time-dependent self-similar 
solution, on the other hand, the fact that timely- dependent 
behavior of the magnetic and gas pressures becomes same is 
one of limits the self-similarity solution. On the other hand, 
the physical quantities with a same physical dimension have 
similar behaviors in self similar solution. 



3.2 Asymptotic behavior 

In this subsection, the asymptotic behavior of the equations 
(22), (26)- (30), and (34) at rj ->• and 7 < 5/3 is investi- 
gated, the asymptotic solutions are given by 



R{rj) ~ RoV~ 3/2 
P(rj) ~ P r,- 5/2 
V(rj) ~ V rj- 1/2 
u)(rj) ~ cjo?T 3/2 
B(rj) ~ B r,- 1/2 
rh(rj) ~ —4ttRoVo 
P{rj) ~ (B 2 /P )rj^ 2 
in which 
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The achieved results for asymptotic behavior of physi- 
cal quantities show that the physical quantities of accretion 
flow are very sensible to parameters of a, 7, /, /3j n , and rhin- 
The Pi„ and rhin are amounts of /3 and rh at rj in that r)i n 
is a point near of the center. The affects of the viscous pa- 
rameter a and the advection parameter / on accretion flow 
are plotted in figure 1. The angular velocity profiles indicate 
that by increasing the viscous parameter a, the angular ve- 
locity of accretion flow decreases, because we increase the 
viscous torque by increasing parameter a. Also increasing 
the advection parameter / decreases the angular velocity 
that is qualitatively consistent with AF06. Figure 1 shows 
the radial infall velocity increases by adding a and / that are 
similar to the results of AF06 and KF08. Also the density 
profiles represent density decreases by adding / and a. 

3.3 Numerical solutions 

If the value of r)i n is guessed, that is a point very near to the 
center, the equations can be integrated from this point to the 
outward through the use of the above expansion. Examples 
of such solutions are presented in figures 2, 3, and 4. The 
profiles in figure 2 are plotted for different p in , the profiles in 
figure 3 are plotted for different / and in figure 4 transonic 
behavior of the accreting gas for different amount of / and 
f5i n is considered. The delineated quantities (Log(-q 3 ^ 2 R), 
Log(—r] 1 / 2 V), ...) in figures 2, 3, and 4 are constant in steady 
self-similar solutions (Narayan & Yi 1994; Narayan & Yi 
1995; Shadmehri 2004; AF06; Ghanbari et. al. 2007; Abbassi 
et al. 2008), while here, they vary by position. 

Figure 2 informs us that density and the radial thickness 
of disk decreases by adding strength of the toroidal magnetic 
field, these results are well consistent with KF08. Also, by 
decreasing amount of magnetic field, the behavior of den- 
sity becomes similar to non- magnetic case (Ogilvie 1999). 
The behavior of the gas pressure in KF08 had polytropic be- 
havior and this selection caused the gas pressure follow the 
density behavior, while here we see behavior of the gas pres- 
sure does not follow the density behavior. Also, by adding 
the P parameter, the radial infall velocity increases; such 
property is qualitatively consistent with AF07 and KF08. 
This is due to the magnetic tension terms, which dominate 
the magnetic pressure term in the radial momentum equa- 
tion that assist the radial infall motion. The profiles of the 
angular velocity imply that the disk is sub-Keplerican in in- 
ner part of the disk and is super-Keplerian in outer part 
of it, while in polytropic accreting flow (KF08) and non- 
magnetic accretion flow (Ogilvie 1999) the angular velocity 
is sub-Keplerian in all radii (KF08). Similar to the results 
KF08 the P parameter, the ratio of the magnetic pressure 
to the gas pressure, is a function of position and arises from 
inner to outer that the result is well consistent with obser- 
vational evidence obtained by some authors (Aitken et al. 
1993; Wright et al. 1993; Greaves et al. 1997). While the j3 
parameter in steady self-similar solution becomes constant 
at all radii (AF06) that is one of restriction of steady self- 
similar solution. Figure 3 is plotted for different amounts of 
the advection parameter /. The advection parameter / has 
slight effect on the toroidal magnetic field, the parameter of 
P, and the Mach number, however has outstanding effect on 
the density, the gas pressure, the radial infall velocity, and 
the angular velocity. The density and the radial thickness 



of disk decrease by more advecting of accreting gas that is 
same at all part of the disk, the result can be achieved by 
assuming of / as a constant amount. Also we see by in- 
creasing the amount of the advection parameter /, the gas 
pressure decreases. By increasing /, the radial infall velocity 
increases and the angular velocity decreases. The results are 
qualitatively consistent with the results of AF06. 

The Mach number profiles in figure 4 imply that the 
flow of outer part for all selected amounts of the magnetic 
field become super sonic. We can see this result in polytropic 
accretion flow by KF08. The advection parameter decreases 
the amount of the Mach number slightly. 

The profiles of physical quantities in figure 2 imply that 
they have the power of law dependency to r\ in magnetical 
domination (p in > 1). So, by fitting a power function on 
data in magnetical domination (p in = 10), we can write 
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The achieved results are different with steady magnetical 
dominated accretion flow (Meier 2005, Shadmehri & Kha- 
jenabi 2005). 



4 SUMMARY AND DISCUSSION 

In the paper, the equations of time-dependent of advection 
dominated accretion flow with a toroidal magnetic field have 
been solved by semi-analytical similarity methods. The flow 
is not able to radiate efficiency, so we substituted the energy 
equation instead of polytropic equation that KF08 had used. 
A solution was found for the case 7 < 5/3 that has differen- 
tial rotation and viscous dissipation. The flow avoids many 
of the strictures of steady self-similar solutions (Narayan & 
Yi 1994; AF06; Ghanbari et al. 2007; Abbassi et al. 2008). 
Thus, the radial-dependence of calculated physical quanti- 
ties in this approach are different from steady self-similar 
solution. 

Increase of the advection parameter / and the param- 
eter Pin will separately increase the infall radial velocity 
and decrease the angular velocity. The flow has differential 
rotation and is sub-Keplerian in inner part and is super- 
Keplerian in large radii in which the behavior is seen in some 
astrophysical objects such as M81, M87 and Milky Way (So- 
fue 1998; Ford & Tsvetanov 1999). The solution showed that 
the flow for all selected amounts of / and Pi„ becomes su- 
per sonic in large radii and sub-sonic in small radii that are 
qualitatively consistent with the results of KF08. The pa- 
rameter of P is a function of position that raises from inner 
to outer and states the magnetic field is more important in 
large radii. It is also consistent with observational evidences 
in the outer regions of YSO discs (Greaves et al. 1997; Aitken 
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et al. 1993; Wright et al. 1993) and in the Galactic center 
(Novak et al. 2003; Chuss et al. 2003). 

Here, latitudinal dependence of physical quantities is 
ignored, while some authors showed that latitudinal depen- 
dence is important in the structure of a disk (Narayan & Yi 
1995; Ghanbari et. al. 2007). Latitudinal behavior of such 
disks can be investigated in other studies. Also we did not 
consider relativity effect, If the central object is relativistic, 
the gravitational field should be changed. Furthermore, in 
a realistic model the advection parameter / is a function of 
position and time, other researchers can consider such disks. 
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Figure 1. Numerical coefficient ujq (dotted lines), Rq (solid lines) and Vb (dashed lines) as functions of advection parameter / or the 
the viscous parameter a. The ratio of specific heats is set to be 7 = 1.5 and the inner mass accretion rate is rhi„ = 0.001. 






Figure 2. Time-dependent self-similar solution for 7 = 1.5, a = 0.5, / = 1.0, and rhi n = 0.001. The lines represent /3; n = 0.1, 0.5, 1.0, 10 
that 0i n is value of /3 in r)i n . 
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Fi gure 3. Time- dependent self-similar solution for '*y — 1.5, en. — 0.5, (3%n — 1.0, and yfii n = 0.001. lines represent / = 0.1, 0.5, 1.0. 




Figure 4. Left panel: Mach number profiles for 7 = 1.5, a = 0.5, / = 1.0, and rhi n = 0.001. Right panel: Mach number profiles for 
7 = 1.5, a = 0.5, f3 ln = 1.0, and m in = 0.001. 



